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Abstract 

The isospin mass splittings of the pseudoscalar and vector D and B light- 
heavy quark system have been calculated using the method of QCD sum rules. 
Nonperturbative QCD effects are shown to be very small, so that mass splittings 
arise almost completely from current quark mass splitting and electromagnetic 
effects, for which a new gauge invariant QED treatment is used. The results 
are consistent with experiment. A measurement of the isospin splitting of the 
vector B mesons would give valuable information about quark mass splittings. 
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1. Introduction 



The QCD sum rule method is well suited for the study of the mass splittings of isospin 
multiplets. In principle, this method allows a consistent treatment using QCD and 
QED field theory of the three sources of isospin symmetry breaking: the current quark 
mass differences in the QCD Lagrangian, the nonperturbative QCD isospin violations 
which arise from u-d flavor dependence of vacuum condensates, and electromagnetic 
effects. In our first work on the isospin splittings of heavy-light mesons]]]], however, we 
found that the standard treatment using the two-loop diagrams for electromagnetic 
corrections to order a e was not gauge invariant. Recently we have developed a gauge- 
invariant theory for QED corrections to the two-point functions used in QCD sum 
rules0, which permits a consistent treatment of all sources of isospin splitting within 
the QCD sum rule method. 

In the present paper we obtain the mass splittings of the D and B pseodoscalar and 
vector mesons. Among the first applications of the method of QCD sum rules was the 
study of isospin violations in the p — uj system (see Ref.||] for a review of the early 
work in this area). Recently, the technique has been used to study the neutron-proton 
mass difference^ |6|, the octet baryon mass splittings |7[] and the mass differences in 
the charmed meson systems (the D and D* scalar and vector mesons) || . In contrast 



with the phenomenological models of heavy-light meson isospin violations || [TO], [TT 
the method of QCD sum rules allows one to separate the current quark mass from 
the nonperturbative QCD effects. Since the nonperturbative QCD isospin violations 
are smaller for the heavy-light quark systems than for the light quark systems [|I], 



12| , the isospin splittings in these charmed and bottom meson systems should be 
given largely by quark mass differences and electromagnetic effects. Moreover, since 
the electromagnetic effects and quark mass splittings enter with different relative 
signs in the charmed vs the bottom systems, the isospin splittings for a heavy-light 
quark system are particularly promising for a determination of the current quark 
masses. (See Ref.Q for a discussion of early work on isospin violations and quark 
mass differences.) 

The formulation of QED corections given in our previous publication ||, in which 
we solved the gauge invariant problem by including an additional current proportional 
to the total charge of the charged current as required by the expansion of the link 
operator, provides a consistent framework to calculate the gauge invariant two-point 
functions in the operator product expansion. Using this theory we find that the 
method is consistent with experiment, which was not true in our earlier attempt ||T| 
to use phenomenological Coulomb potentials to estimate electromagnetic effects; and 
that we can obtain a consistent QCD/QED sum rule analysis of the isospin mass 
splittings of the D,D*, B and B* systems. This will be our main conclusion. 

The starting point of the QCD sum rules is to evaluate the two point function 



IV(g 2 ) = z J d 4 xe^ < T(J^x)J u (0)) > 



2 



in the Wilson operator production expansion. To investigate the violation of the 
isospin symmetry, the two point function H(q 2 ) can be written as 

n(g 2 ) = n (g 2 ) + m q U m (q 2 ) + U em (q 2 ). (2) 

Uo(q 2 ) in Eq. is the leading term for the light heavy quark system, the isospin viola- 
tion comes from the isospin splitting of the quark condensates in the non-perturbative 
power correction, which has been given in previous publications]!], The second 
term in Eq. |2] is the light quark mass, m q , expansion of the total correlator; and 
the third term corresponds to the electromagnetic corrections, which have been not 
treated consistently in the literature. 

The paper is organized as follows. The light quark expansion correlator, U m (q 2 ), 
and the electromagnetic effects, IT em , are discussed in Sections 2 and 3. In Section 
4, we shall present our numerical analysis of the isospin splittings for the charm and 
bottom mesons, and finally the conclusion is given in Section 5. 



2. The Contributions From The Up and Down Quark Mass 
Difference 

We write U m (q 2 ) in Eq. ^| as 

U m (q 2 ) = Cjl + C 3 (qq) + C 5 (q(° " G)q), (3) 

where the coefficient Cj represents the perturbative contributions and we limit the 
expansion to dimension = 5 operators. Based on the results in Ref. |14[], we find the 
m q expansion to order ot s is 
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and l(x) — — Jq ln(l —y)^ is the Spencer function. 

The coefficients C3 and C5 in Eq. |^ come from expansions in the small quark mass 
of the light quark propagator. The relevant terms are 
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where x = 7 M x M . The result to the order a s for the coefficient C3 is 
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for the vector currents. The m q expansions of the one loop corrections to the quark 
condensate are finite after the mass renormalization, in which the pole mass is used. 
The coefficient C5 in Eq. 2 is 
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for the pseudoscalar and vector currents, respectively. 



(10) 



3. The Electromagnetic Effects 

We write the electromagnetic corrections in analogue to Eq. [3]; 

n em (g 2 ) = DjI + D 3 (qq) + D 5 (q(° ' F)q), (12) 



where the first term is a two-photon-loop perturbative correlator, the second term 
represents a photon-loop correction to the quark condensate, and the third is a di- 
mension 5 quark-photon condensate in analogue to the quark gluon condensate in 
QCD. In reality, the quark-photon condensate, (q(a ■ F)q), is not well understood 
because of the gauge invariance problem, although it has been used in the analysis 
of the proton and neutron mass difference 0. Thus, we introduce a new parameter 5 
by relating the quark-photon condensate to the quark-gluon condensate; 

(q(a-F)q) = -5(q(a-G)q), (13) 

where F and G correspond to photon and gluon fields respectively. 
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The gauge invariant electromagnetic current at order a e is 

J^(x) = q{x)lv,Q{ x ) +«e T g(x)7 M / A u (y)dy v Q(x) 1 (14) 



o 

where the second term is a gauge fixing term for the charged currents, and comes 
from the expansion of the link operator exp(ieT f% A v (z)dz 1/ ) which is being inserted 
between the heavy and light quark fields to make the operator product expansion in 
QED gauge invariant, and is the total charge of the system. 

Thus, the operator product expansion in QED for charge neutral currents at order 
a e should have the same form as that in QCD at order a s except that the coupling 

constant 4^- in QCD is replaced by - 2 — - in QED (note that for neutral currents the 
charges for the light and heavy quarks are equal with opposite sign). By separating 
the gauge dependent and independent terms for the leading current J M = q(x) r Y li Q(x), 
we have shown in Ref. that it is convenient to write the operator product expansion 
in QED into three parts; 

D C = D° + D Q + D q . (15) 

The first term, D°, in Eq. [15] is the gauge independent term for the leading current 
J ^ = q(x)jnQ(x), thus, it also has the same form as its counterpart in QCD except 
that the coupling constant is replaced by — £Qe ^ ae in this case. The second and 
third terms in Eq. [15| include the gauge dependent terms of the leading current that 
only involves the self-energy diagrams of the heavy and light quarks and the terms 
generated by the gauge fixing current, J'^ = iexq^x)^^ Jq A u (y)dy u Q(x), in Eq. [14]. 

The advantage of this approach is that the first term in Eq. |l| and the correlators 
for the charge neutral currents have been derived in QCD[[L4]]; we have 
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for the charge neutral pseudoscalar and vector currents respectively, where f{x) is 
given in Eq. |5]. 



For the charged current, the first term in Eq. 03 is essentially given by Eqs. 16 



and [T7| for pseudoscalar and vector mesons, and we found that D q in Eq. [T5| is zero 
since it is proportional to the wavefunction renormalization for a zero mass particle in 
the Landau gauge 0. Therefore, the only term left in Eq. [1^ is D®, and its imaginary 
part is shown to be0 
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where = e g + eg for charged systems. 

Similarly, one can calculate the one-photon loop corrections to the quark conden- 
sate. The coefficients D3 for the charge neutral system are just an extension of the 
one loop corrections in QCD; therefore, we have 



e Q a e 



M, 



Q 



2tt q 2 (l-x) 



1 + 



-3(1 -a;) (l+zln 



1 — x 



(20) 



and 



e 2 Q a e 



Mr 



2tt q 2 {l 



x 




x 



x In 



x 



x 



(21) 



for the charged neutral pseudoscalar and vector currents respectively. 

The coefficient D3 for the charged currents follows from the same procedure used 
in deriving the coefficient Dj. The first term in Eq. [15], D® in this case, can be 
obtained from Eqs. ^0] and ^1] by changing the coupling constant to — tQt q a e . We 
find that becomes zero after the mass renormalization, in which the pole mass 
is used. This is because Df only corresponds to the self-energy diagram, and its 
wavefunction renormalization integral is canceled exactly by the gauge fixing term, 
for which the physics is the same as that of the vanishing D q for the two-photon-loop 
diagram. The term is found to be 
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(22) 



for both pseudoscalar and vector currents. It is interesting to note that the term D3, 
present in both charged and neutral mesons, has opposite signs for the pseudoscalar 
and vector currents. It plays the same role as the hyperfine interaction in the phe- 
nomeno logical quark models, while _D§, present only in the charged current, could be 
interpreted as the contribution from the photon cloud. This shows the subtle corre- 
lation between the photon interaction and the light quark condensates. In particular, 
the effects for the charged bottom mesons should be larger than the charged charm 
mesons, as the factor eTe q for B mesons in Eq. |22] is two times larger than that of 
the D mesons. 
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Similar to the quark-gluon condensate, one could obtain the quark-photon con- 
densate term in the fixed point gauge, and there is no contribution from the gauge 
fixing term in Eq. |TJ| to the coefficient D 5 . Therefore, the coefficient D 5 is a simple 
extension of the corresponding QCD diagrams, and we have 

5 " V e W 2(g 2 - Mq) 3 + ^ e Q e q ) - Mq) 2 ^ 

and 

where the top and bottom matrix elements correspond to the neutral and charged 
currents respectively. Eq. [24] and the first term in Eq. ^ come from the term in 
which a photon is being emitted by the light quark itself, whose counterpart in QCD 
is the fifth term in Eq. [7| They have the same effects as the light quark mass differ- 
ence in B and D mesons, and physically could be regarded as the contributions from 
the electromagnetic effects to the light quark masses. The second term in Eq. [23] 
comes from the photon being emitted from the heavy quark, and thus has different 
signs for the charged and neutral systems. The physics of this term is the same as 
the hyperfine interaction in the phenomenological quark models, which generates the 
mass differences between the vector and pseudoscalar states. Therefore, the domi- 
nant electromagnetic effects should come from the photon corrections to the quark 
condensate and the two-photon loop corrections, while the quark-photon condensate 
term is only responsible for the difference between the vector and pseudoscalar states. 

4. The Numerical Results 

Defining the quantity uj 2 = Mq — q 2 for the Borel transformation, where uj 2 measures 
the off-shell effects of the heavy quarks, we find 

4 Mp-MQ 

fp^'^ = (25) 



and 



= W ^b) (26) 

Q 

for the pseudoscalar and vector current. The correlator U t (u^) can be divided into 
the perturbative and nonperturbative parts: 

w 2 

U t (u%)= f ImCj(uo 2 )e~^duj 2 + Tl np (uj 2 B ), (27) 
Jo 

where the perturbative term is written in dispersion integral form, and the nonper- 
turbative correlator n np (a;^) comes from the standard Borel transformations of the 
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terms associated with the condensates. The masses M p and M v in Eqs. |12] and [16 
can be obtained from 

M^+ ^y . (28) 

The mass differences between different isospin states are obtained by adjusting the 
parameters s° in Eq. [27] to insure that the mass difference is independent of Borel 
parameter u B . The range of io 2 B is chosen in such way so that the continuum and 
the nonperturbative contributions are minimum. By fitting the masses and decay 
constants we find the working range for the Borel parameter u 2 B are 1.4 < uj 2 b < 2.4 
for D mesons, 0.9 < uo 2 B < 1.9 for D* mesons, 3.9 < uj 2 b < 4.9 for B mesons, and 
3.3 < u B < 4.3 for B* mesons, where u B has a unit of GeV 2 . 

In order to fit both D and B systems in the same time, we also include the 
dependence of the light quark mass and the quark condensate on the factor, 

L a = [ln(M Q /A)/ln(/i/A)r (29) 

where a = — ~ for light quark masses and | for the quark condensate are the anoma- 
lous dimensions. The quark condensate (qq) has the standard value — 0.0138GeV 3 , 
and the quark-gluon condensate is related to the quark condensate by (qa ■ Gq) = 
ml(qq) with m 2 , = 0.8 GeV 2 . As a first approximation, we let m\ be independent 
of the isospin splitting so that the ratio of quark-gluon condensates between up and 
down quarks is the same as that of the quark condensate, which is determined by the 
parameter 7; 

7=^-1- (30) 
(mm) 

There are some constraints on the isospin differences of the light quark masses 
and condensates: the estimate from the chiral perturbation theory gives ma — m u = 
3.0 MeVQ, while the parameter 7 in various estimates ranges from —0.00790 to 
—0.0020, ||] . However, the parameter 7 in previous studies may not be reliable 
because the electromagnetic effects have not been consistently included. As a first 
attempt, we set md — m u = 4 MeV, 7 = —0.002, and 5 = 0, by adjusting the threshold 
parameter s° in Eq. ^ to obtain a constant in the working region of lj b , the resulting 
isospin mass splittings are 



and 



in units of MeV. The contributions from the electromagnetic effects to the isospin mass 
splittings can be obtained from the difference between the isospin mass splittings with 
and without the electromagnetic effects, and results are approximately 

\m ± -M°\ =0.4 (33) 

I- J em 



s 



for both D and D* mesons and 



{ M ± — M°\ =1.9 (34) 

l J em 

for both B and B* mesons. Thus, the contributions from the electromagnetic effects 
are considerable larger in B mesons. The dominant terms in the electromagnetic 
contributions are the photon loop corrections to the quark condensate, and the quark 
condensate for the B mesons is also considerable larger than that for D mesons due 
to the factor L a in Eq. |29| in addition to the electric charge factor in Eq. |22|. 

It is interesting to consider the phenomenological values for the electromagnetic 
corrections obtained in previous work on isospin splitting of the heavy-light mesons. 
In the quark model estimate of Ref. in which the {M ± — M°} em is taken as the 
Coulomb potential, the values for the D and B systems are approximately equal. 
This is not satisfactory, since {M ± — M°} em must be much larger for the B systems 
than the D systems to agree with experiment. In the work of Ref.|TT| {M^ — M°} em 
for the D systems is found to be 7/5 times that of the B, causing an even larger 
discrepancy with experiment if used by us. 

We find little dependence of the isospin mass splitting on the parameter 7; varying 
the parameter 7 from —0.002 to —0.004 only generates a 0.1 MeV increase of the mass 
splitting for D mesons, and no significant increase at all for B mesons. However, this 
may be due to the assumption that the ratio of the quark gluon condensates has 
the same dependence on the parameter 7 as that of the quark condensates between 
the up and down quarks. This shows that the contributions from the quark gluon 
condensate is important, and the heavy-light quark systems may not be a good place 
to determined the parameter 7. 

We show our best fit to the isospin mass splittings in Fig. 1 for D, Fig. 2 for D*, 
Fig. 3 for B, and Fig. 4 for B* states. The parameters for this fit are — m u = 3.0 
MeV, 5 = 0.01, and 7 = —0.002. The resulting isospin splittings are 



and 



The corresponding differences of the threshold parameters s are 



and 



± o _ 0.012 for D 

I 0.026 forD* 1 > 



s± s °-\ _o.Q26 forfi* (38) 
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in units of GeV 2 , which indeed have the same signs as those of the light quark mass 



difference ||]. This is in good agreement with data (M D ± — M D o = 4.7 ± 0.7 MeV[16] 



and M D *± - M D *o = 3.32 ± 0.13 MeVM for D mesons and -0.1 ± 0.8 MeVfliq for 



B mesons). Perhaps more encouraging is that the light quark mass difference in this 
fit is in good agreement with the result of the chiral perturbation theory ||. 

An important feature of this investigation is the differential vector-pseudoscalar 
isospin mass splittings A(M); 

A(M) = (M* ± - M*°) - (M ± - M°), (39) 

which provides very important probe of the hyperfine splittings of the heavy-light 
quark potential. The quantity A(M) comes from two sources: the hyperfine inter- 
action from the differences of the light quark masses and condensates and from the 
electromagnetic effects (in particular, the quark-photon condensate in Eq. f2^). These 
two sources have the same sign for the D mesons and opposite signs for B mesons. 
Therefore, the quantity A(M) for B mesons really depends on the competition be- 
tween these two sources. From Eqs. |35| and [36], we have 



A/ ,,x f —0.8 for D mesons 

AM = <^ 40 
0.1 tor B mesons 



This is in good agreement with recent studies by Cutkosky and Geiger[[0|, the quan- 
tity A(M) for B mesons is essential zero within the theoretical uncertainty. A mea- 
surement of the isospin mass splittings for the vector B mesons will provide a crucial 
test in this regard. A study in chiral perturbation theory finds ]ll| that the contri- 



butions from quark mass difference only gives 0.3 MeV for D mesons, and the rest 
might comes from the electromagnetic corrections. Our calculation is consistent with 
this conclusion, and furthermore we also find that there is only a small contribution 
from the difference of the quark condensates to A(M). 



5. Summary 

Gauge invariance dictates that there is an additional gauge fixing current for charged 
mesons. This gauge fixing current generates additional contributions to the electro- 
magnetic effects for charged mesons. In particular, we find the photon-loop correction 
to the quark condensate dominant for electromagnetic effects, and because of this cor- 
relation between quark condensates and photons, the electromagnetic corrections are 
considerably larger in B mesons than those in D mesons. This is quite contrary to 
the conclusions of phenomenological quark models. Indeed, our results are in good 
agreement with the known data, and further experimental studies of B* mesons will 
certainly provide us more insights into the hyperfine interactions in the quark-quark 
potential. 

We find little dependence of the isospin mass splitting on the isospin differences 
of the quark condensates, a nonperturbative QCD effect given by the parameter 7 in 
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our approach, while the light quark mass difference used in our calculation is in good 
agreement with the result of the chiral perturbation theory. Therefore, a consistent 
treatment of the electromagnetic effects would allow us to study the light quark mass 
difference in heavy-light quark systems. 

Another interesting application of our approach is the isospin splittings in Kaon 
systems. Since the strange quark has much less mass than the charm and bottom 
quarks, the mechanism of the electromagnetic effects might be different from the 
heavy light quark systems. This investigation is in progress, and the result will be 
published elsewhere. 

This work is supported by National Science Foundation grant PHY-93 19641. 
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Figure Captions 

1. The sum rule for D mass splittings: M^ — M , with parameters m^ — m u = 
3.0 MeV, 5 = 0.01 and 7 = -0.002. The difference of the threshold 
parameter is s — s° = 0.012 GeV 2 for the center curve, 0.016 and 0.008 
GeV 2 for the upper and lower curves repectively. 

2. The same as Fig. 1 for D* mass splittings with the threshold parameter 
difference 0.026 for the center curve, 0.03 and 0.022 GeV 2 for the upper 
and lower curves respectively. 

3. The same as Fig. 1 for B mass splittings with the threshold parameter 
difference 0.024 for the center curve, 0.032 and 0.016 GeV 2 for the upper 
and lower curves respectively. 

4. The same as Fig. 1 for B* mass splittings with the threshold parameter 
difference 0.026 for the center curve, 0.034 and 0.018 GeV 2 for the upper 
and lower curves respectively. 
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